Theory of the hourglass dispersion of magnetic excitations in high-T c cuprates 
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A theory for the dispersion of collective magnetic excitations in superconducting cuprates is 
presented with the aim to cover both high and low doping regimes. Besides of spin fluctuations 
describable in the random phase approximation (RPA) we allow for local spin rotations within a 
mode-coupling theory. At low temperatures and moderately large correlation lengths we obtain 
two branches of excitations which disperse up- and downwards exhibiting the hourglass behavior 
observed experimentally at intermediate dopings. At large and small dopings our theory essentially 
reduces to the RPA and spin wave theory, respectively. 
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The low-temperature magnetic response of many high- 
Tc superconductors is characterized by a resonant mode 
inside the superconducting gap around the antiferromag- 
netic wave vector Q. This collective mode manifests it- 
self as a single peak at Q which splits into two peaks 
dispersing up- and downwards in frequency away from 
Q. This unusual dispersion resembles the shape of a 
hourglass [H 0]. Theories to explain thisphenomena use 
either a more localQ or an itinerant [3, [f| description of 
the magnetism. The second approach considers particle- 
hole excitations with spin flips which interact within the 
random phase approximation (RPA) forming a dispersing 
bound state in the superconducting gap. This approach 
yields only one branch of excitations below the Stoner 
continuum whereas it has been established recently that 
the lower branch, the center of the hourglass as well as 
part of the upper branch lie below this continuum in the 
gapped region A more theoretical argument for the 
incompleteness of a RPA description comes from the fact 
that different spin directions do not mix as a function 
of time in this approximation which excludes local rota- 
tions of spins known from spin wave theory. Below we 
will present a theory which contains both spin wave the- 
ory and RPA as special cases. At intermediate dopings 
we will show that both RPA and spin wave like spin fluc- 
tuations are important and produce the two branches of 
the hourglass dispersion. 

We consider the t-J modelQ with the Hamiltonian H, 

ff = ^e(k)cLc kCT + ^J(k)S k S_ k . (1) 

kcr k 

c ko .,6k CT ar e creation and annihilation operators, respec- 
tively, for electrons with momentum k and spin projec- 
tion a excluding any double occupancies of sites. S k are 
spin operators in momentum space, e(k) and J(k) are 
the bare electron dispersion and the Fourier transform 
of the Heisenberg coupling, respectively. A useful ap- 
proximation for H, used in the following, is obtained by 
taking the large N limit of the t-J model, where e(k) de- 
scribes a renormalized dispersion of quasi-particles and 



the fermionic operators can be treated as usual creation 
and annihilation operators. 

In the following we are interested in the time evolution 
of the spin operator S k = 1/2 A(k, q), 

A(k, q) = Ck+qa^/JCq/^ ( 2 ) 

a, 13 

where <r denotes the vector of the three Pauli matrices. 
It obeys the equation of motion 

^ = <(A) + A)A, (3) 

with 

£ A(k,q) = (e(k + q)-e(q))A(k,q) 

+J(k)(/(q)-/(k + q))S kl (4) 



AA(k, q ) = -i^J(kO- 
k ' 

(S k , X (A(k - k', q) + A(k - k', k' + q))) . (5) 

/(k) is equal to (c k(7 c kcr ), where (...) denotes the thermo- 
dynamic expectation value, and x stands for the vector 
product. Since we are only interested in the spin re- 
sponse we have dropped terms on the right-hand side 
of Eq.© which involve fluctuations in the density. We 
also dropped an overall prefactor denoting the number of 
primitive cells. The unperturbed Liouville operator Cq 
describes the time evolution of the system in the RPA. 
From its explicit expression in Eq. ((4|) follows that it does 
not mix different cartesian components of the spin opera- 
tors. In contrast to that the time evolution described by 
C\ involves product states of spin operators, mixes dif- 
ferent spin components and thus can describe rotations 
of spins due to fluctuating fields. 

The spin susceptibility x(k, z) can conveniently be cal- 
culated from the associated Kubo relaxation function 
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$(k, z) = (x(k) — x(k,^))/z where 2; is a complex fre- 
quency and x(k) is equal to x(k, 2 = 0). Due to the 
rotational invariance in spin space we may assume that 
X, M etc. always refer to the z-direction. Using the Mori 
formalism $ can be written asQ 



$(k,z) 



X(k) 



2 + M RPA (k, z) + M(k, z) 



(6) 



The first memory kernel M RPA (k, z) describes the time 
evolution of spin operators by C - According to Eq.([31) 
the direction of the spin operators is conserved and they 
remain always linear in the operators A. Eliminating 
the A operators in the equation of motion in favor of 
the orginal S operators yields an explicit expression for 
M RPA which may be expressed in terms of the RPA spin 
susceptibility x RPA (k, z), 



M 



with 



RPA/ 



(k, z) = z X KPA (k, z)/( X « PA (k)- X ^ A (k, z)), (7) 



l/ X RPA (k, W ) = l/x (0) (k^) + J(k), 



(8) 



X^(k, u>) is the free susceptibility. 

The second memory contribution M (k, z) is due to 
the time evolution of single spin operators S into prod- 
uct states of spin operators described by C\. Using the 
mode-coupling assumption and performing the analytic 
continuation z — ► uj + ir), we obtain for the imaginary 
part of M, denotd by M", 

M"(k,w) = ^(J(k')- J(k-k')) 2 ^"(k,k', W )/(o;x(k)), 

k' 

(9) 

£>"(k,k',w) =tt J duj'A(k-k',u-Lj')A(k',uj') 

(b(Lu')-b(uj' -LU)). (10) 

A is the spectral function of the spin propagator and b(u>) 
the Bose function. 

SMALL CORRELATION LENGTHS 

In the case of small antiferromagnetic correlation 
lengths £, corresponding to the overdoped regime, the 
RPA should be a reasonable approximation for the spin 
susceptibility. The dashed line in Fig. 1 shows the imagi- 
nary part of X RPA (Qi w ) f° r Q = (^i 7r ) using the param- 
eters tb2 in Table I of Ref.[4( and a chemical potential 
corresponding to the doping 6 = 0.20. The energy unit is 
1 eV in the following and lengths are measured in units of 
the lattice constant a of the square lattice. The supercon- 
ducting order parameter is A(k) — A(cosfc :E — cosfc y )/2 
with A = 0.029, J equal to 0.135 and r] = 0.004. The 
dashed line in Fig. 1 illustrates that most of the spectral 




FIG. 1: x (Q> w ) without (dashed line) and with (solid line) 
memory function M for a doping 8 = 0.20 and T = 0. The 
dashed curve corresponds to the RPA. Inset: real part M' 
and imaginary part M" of M(Q,u;). 

weight resides in the bound state at the energy 0.038 and 
that only a small part of it has been left in the contin- 
uum at higher energies. Away from (jr, n) the dashed 
curve in Fig. 1 does not change dramatically as long as 
the bound state lies still in the gapped region. Entering 
the particle-hole continuum by going further away from 
(tt, tt) destroys the bound state and most of the spectral 
weight shifts to high energies of the order of t. The pa- 
rameters used in Fig. 1 yield £ ~ 0.78. For such a small 
£ practically all momenta in the sum over momenta in 
Eq.© contribute substantially which means that M" is 
mainly determined by contributions away from the small 
region around (tt, tt) so that the bound state and its low- 
energy scale is rather irrelevant for M. This is confirmed 
by an explicit calculation of M using RPA results for the 
various quantities in Eq.@. The result is shown in the 
inset of Fig. 1 for T=0. M" (solid line) is structureless 
except at small energies where it vanishes rapidly due to 
the smallness of A in this region and the cutoff u> for the 
integration over ui' in Eq. (fTU|) due to the bosonic factors. 
Taking M into account in Eq.® yields the solid line in 
Fig. 1 which differs only marginally from the dashed line. 
This shows that at short correlation lengths the RPA re- 
sult for x" is essentially correct and that the correction 
M to M RPA is rather small. The underlying physical pic- 
ture is that the momentary local axis of preferred spin di- 
rections fluctuates very rapidly due to the random forces 
induced by C\. The spectrum of these forces is given 
by M" and characterized by the large energy scale t in 
agreement with the inset of Fig. 1. 



LARGE CORRELATION LENGTHS 

For large £ the spectral function A(\s.',u>') is strongly 
peaked at k' = Q. This means that the integration 
over k' in Eq. ^ is restricted to momenta near or near 
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Q. Since we are interested in external momenta k ~ Q 
the momentum of one of the two spectral functions in 
Eq. tfTT))) is small. Due to spin conservation this spectral 
function describes spin diffusion and is mainly restricted 
to small values of a/. As a result one may neglect the 
small frequency transfer in the second spectral function 
in Eq. (fT0|) . Taking also the real part of M into account 
we obtain from Eqs.© and (fT0|) . 



M(k,^) = -^ 2 (k)<E>(Q,^)/x(Q), 



(11) 



with 



^ M = ~k\ £ (J(q) - J ( k - q)) 2 (^- q Vk)x(q)), 

X( k ) q 

and the equal-time correlation function 



(13) 



(SkS-k) = / du;b(u;)A(k,u;). 



In deriving Eq. ljlip we used the fact that the two memory 
functions in Eq. ([6]) depend for our parameters only slowly 
on momentum around the wave vector Q so that the com- 
bination <I>/x on the right-hand side of Eq. fTTj) may be 
evaluated at Q. The sum over q in Eq.(fT2"j) runs over 
half of the Brillouin zone centered around Q. The evalu- 
ation of the above expressions using the RPA encounters 
a problem: £, calculated in the RPA, is in the optimal 
and moderately underdoped region around one or smaller 
and increases substantially only near the transition to the 
antifcrromagnetic state in disagreement with the experi- 
ment. For instance, we have for 5 — 0.12 £ ~ 0.8, Ref.[5| 
£ ~ 0.6 using quite different parameter values, whereas 
the experimental values for £ are larger by about a fac- 
tor 5 or more[8j. Since this large discrepancy would af- 
fect severely the momentum sum in Eq. (1121) we prefer to 
use a realistic xO*) as input in calculating M and write 
X(k) = x(Q)/(l+£ 2 (Ak) 2 ) for Ak = k-Q - consider- 
ing £ as a parameter to be determined from experiment. 

It is instructive to study the frequency dependence 
of the denominator of Eq.©. In order to describe a 
slightly underdoped case we choose the same parame- 
ters as in Fig. 1, a chemical potential corresponding 
to S = 0.12, £ = 5, and the cutoff l/£ for the sum 
over q in Eq. (fT2")) . The solid and dotted line in Fig. 
2 show 1 + Af' RPA (k,w)/w for H=0.5 and 0.38, respec- 
tively, writing k = (H, 0.5)27r. This quantity is prac- 
tically independent of momentum, increases monotoni- 
cally with u> and is zero at the RPA resonance energy 
uir ~ 0.038. The dashed and dash-dotted lines in Fig. 2 
show — M'(k, uS)/u) for the same momenta. These curves 
resemble the real part of an oscillator located at u>r with 
an oscillator strength being very small at Ak = k— Q = 
and strongly increasing with |Ak|. The poles of Eq.f6|) 
are given by the common points of the two curves de- 
noted by squares and circles. Since the common point at 



k=(H,0.5)2it 

M'(k,a>)to i =0004 
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FIG. 2: (color online) 1 + M' RPA (Q,w)/w (solid and dotted 
lines) and — M'(Q, u)/u) (dash-dotted and dashed lines) as a 
function of frequency for two momenta H. Squares and circles 
denote poles of x(k,oj). 




FIG. 3: (color online) Position of the poles of p<;(k, a>) as 
a function of momentum. The shadowed region marks the 
particle-hole continuum. 

lu = ujji (not shown in Fig. 2) has vanishing pole strength 
there are two branches of collective spin excitations. For 
vanishing damping r\ their dispersion is shown in Fig. 3 
by solid lines. They approximately touch each other at 
Ak = and disperse up- and downwards with increasing 
|Ak|. For not too large |Ak| both branches lie below the 
continuum in agreement with experiment 1]. Performing 
the calculation in the normal state at T = T c the solid 
and dotted lines in Fig. 2 lie everywhere above zero but 
the solid and dashed and also the dotted and dash-dotted 
lines have still one common point at larger frequencies. In 
this case only the upper but not the lower branch exists 
in agreement with experiment 0. At very low dopings 
M RPA — > due to the constraint and the pole condition 
u) + M'(k, u>) yields in the presence of long-range order 
the correct spin wave dispersion^. 

Several prerequisites are necessary to obtain the above 
hourglass dispersion for spin excitations. There must ex- 
ist two different kinds of spin excitations to account for 
the two branches. The first one are RPA spin fluctu- 
ations where all induced spin moments have the same 
direction. They may be characterized by the fact that 
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FIG. 4: (color online) Imaginary part of the spin susceptibility 
at zero temperature as a function of momentum for different 
frequencies. 

the internal fields induced by the Heisenberg interaction 
conserves frequency, momentum and spin direction which 
is a direct consequence of the one-mode behavior of Co in 
Eq.((4]). The second one are local rotations of spins under 
the influence of C\ in Eq.(|5|). In this case a spin in z 
direction acquires in its time evolution also a component 
in x direction due to the presence of a spin fluctuation 
in y direction. The pure form of the two kind of spin 
excitations are obtained for M = and M RPA = 0, re- 
spectively, and are realized approximately at large and 
small dopings. In the hourglass regime M and Af RPA 
are of similar magnitude. 

Another prerequisite for hourglass behavior is that £ 
is substantially larger than 1. Only then is the momen- 
tum integration in Eq. ([9]) restricted to the resonance and 
the spin diffusion modes yielding oscillator-like behavior 
of M near u>r. The up- and downwards dispersion and 
their approximate degeneracy at Q is mainly determined 
by w 2 (k), which according to Eq. (|12p is roughly propor- 
tional to (S , Ak5'-Ak}(£" 2 + (Ak) 2 ). The first factor tends 
to zero at low temperatures for Ak — > and saturates 
at large Ak. As a result w 2 (k) is very small at Ak = 
causing the approximate touching of the two branches 
at oj = lor and Ak = 0. With increasing Ak Lo 2 {k) in- 
creases strongly leading to a downward dispersion of the 
lower branch even if cor was practically dispersionless as 
in our case. According to Fig. 3 the upper branch in- 
creases at large |Ak| roughly as 0.6J|Ak|, i.e., with an 
effective spin wave velocity which is reduced by about 
a factor 2-3 compared to spin wave theory similar as in 
experiment [J Q • 

Using the same parameters as in Figs. 2 and 3 Fig. 4 
shows x" as a function of H with the frequency as a pa- 
rameter. As suggested by Fig. 3 x" exhibits a hourglass 
dispersion with intensities which are largest near u>r and 
decay rather fast and slow towards lower and higher fre- 
quencies, respectively. Fig. 5 shows x" as a function of 
lu for a fixed H as a parameter. In agreement with Fig. 
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FIG. 5: (color online) Imaginary part of the spin susceptibility 
at zero temperature as a function of frequency for different 
momenta. 

4 the strong peak at H = 0.5 splits into two peaks with 
decreasing H which disperse up- and downwards in fre- 
quency. The curve for H = 0.5 calculated for the small 
damping r\ = 0.004 has the shape of a Lorentzian. At 
smaller dampings this peak splits into a double peak due 
to the small gap between upper and lower branch shown 
in Fig. 3. For M = only the lower, weaker peak is 
obtained. 

In conclusion, we have shown that the memory func- 
tion of the spin susceptibility contains in general two dis- 
tinct contributions due to RPA-like and due to rotational 
spin fluctuations. The first one dominates at large, the 
second one at small dopings. At intermediate dopings 
both are of similar magnitude which leads to one upwards 
and one downwards dispersing branch of excitations. At 
low temperatures the two branches are approximately de- 
generate at Q which explains, at least qualitatively, the 
observed hourglass dispersion at intermediate dopings. 

The author is grateful to V. Hinkov, P. Horsch, D. 
Manske and H. Yamase for useful discussions. 



[1] B. Fauque et al, Phys. Rev. B 76, 214512 (2007), and 

references therein. 
[2] V. Hinkov et al., nature physics 3, 780 (2007). 
[3] I. Sega, P. Prelovsek, and J. Bonca, Phys. Rev. B 68, 

54524 (2003). 
[4] MR. Norman, Phys. Rev. B 63, 92509 (2001). 
[5] H. Yamase and W. Metzner, Phys. Rev. B 73, 214517 

(2006), and private communication. 
[6] For a review, see M. Ogata and H. Fukuyama, Rep. Prog. 

Phys. 71, 036501 (2008). 
[7] D. Forster, "Hydrodynamic Fluctuations, Broken Sym- 
metry, and Correlation Functions", W.A. Benjamin, Inc., 

Reading, Massachusetts (1975), chapter 9. 
[8] P. Dai, H.A. Mook, R.D. Hunt, and F. Dogan, Phys. Rev. 

B 63, 054525 (2001). 



